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Abstract 



We report the first observation of color-suppressed B — ► D w , D it , D rj, 
and D°lo decays and evidence for B° — > D*°t] and D*°uj. The branching 
fractions are found to be B(B° -► D°ir°) = (3.1 ± 0.4 ± 0.5) x 10" 4 , B(B° -► 
D* 7T ) = (2.7 ±g;f ±g;|) x 10~ 4 , fi(5° D ^) = (1.4 t° i ± 0.3) x 10~ 4 , 
6(5° -► £>°u;) = (1.8 ± 0.5 ±[J; 4 ) x 10~ 4 , and we set 90% confidence level 
upper limits of B(B° -► D* ??) < 4.6 x 10~ 4 and 5(5° -► D*°u;) < 7.9 x 10~ 4 . 
The analysis is based on a data sample of 21.3 fb" 1 collected at the T(4S) 
resonance by the Belle detector at the KEKB e + e~ collider. 



PACS numbers: 13.25.Hw, 14.40.Nd 



Decay modes such as B° — > D^ + n^ were amongst the first B meson decays to be 
fully reconstructed. However, the B° — > D^*>°h° decay modes, where h° represents a light 
neutral meson, have not been observed to date. A recent search by the CLEO Collaboration 
yielded only upper limits. As shown in Fig. [I], these decays are expected to proceed 
via an internal spectator diagram and be suppressed relative to decays proceeding via the 
external spectator diagram, since the color of the u antiquark produced by the weak current 
must complement the color of the c quark. The contribution of the VF-exchange diagram 
is usually assumed to be negligible @. Studies of such color-suppressed decay modes are 
useful for testing models of hadronic B meson decay and provide information on final- 
state interactions. Final state interactions, if also present in B decays to charmless modes, 
could have significant impact on CP violating rate asymmetries ||. Furthermore, studies of 
B° — > D^°h° where the D° decays to a CP eigenstate will provide access to the CP violation 
parameter sin 20i in b — > cud processes [|[]. Thus, measuring the strength of color-suppressed 
modes can have profound implications on B physics and CP violation. 

In this paper, we report on a search for the color-suppressed B° — > D°h° and D*°h° || 
decays, where the neutral meson h° is either a 7r°, 77, or u. We make first observations of the 
modes B° — > D°ti , D*°tt°, D°t), and D°u and find evidence for the modes B° — > D*°r] and 
D*°uj. The data used in this analysis were collected with the Belle detector [[]] at KEKB J?|, 
a double storage ring that collides 8 GeV electrons with 3.5 GeV positrons with a 22 mrad 
crossing angle. The data sample corresponds to an integrated luminosity of 21.3 fb -1 at the 
T(4f>) resonance, which contains 23.1 million BB pairs, and 2.3 fb -1 taken 60 MeV below 
the resonance. 

Belle is a general-purpose detector with a 1.5 T superconducting solenoid magnet. 
Charged particle tracking, covering 92% of the total center-of-mass (CM) solid angle, is 
provided by the Silicon Vertex Detector (SVD) consisting of three concentric layers of double- 
sided silicon strip detectors and a 50-layer Central Drift Chamber (CDC). Charged hadrons 
are distinguished by combining the responses from an array of Silica Aerogel Cerenkov Coun- 
ters (ACC), a Time of Flight Counter system (TOF), and dE/dx measurements in the CDC. 
The combined response provides K/tt separation of at least 2.5a for laboratory momenta up 
to 3.5 GeV/c. Photons and electrons are detected in an array of 8736 CsI(Tl) crystals (ECL) 
located inside the magnetic field and covering the entire solid angle of the charged particle 
tracking system. The 1.5 T magnetic field is returned via an iron yoke, instrumented to 
detect muons and Kl mesons (KLM). The KLM consists of alternating layers of resistive 
plate chambers and 4.7 cm thick steel plates. 

We reconstruct light neutral mesons using the 7r° — > 77, 77 — > 77, 77 — > 7r + 7r 7r° and 
oj — > 7T + 7r~7r° decay channels. Charged tracks are required to have impact parameters that 
are within ±5 cm of the interaction point along the positron beam axis and 1 cm in the 
transverse plane. We reject tracks that are consistent with being electrons or muons. The 
remaining tracks are identified as pions or kaons according to a kaon to pion likelihood ratio. 
Candidate 7r° mesons are reconstructed from pairs of photons in the ECL if the 7 pairs have 
invariant mass inside a ±3a (a = 5.4 MeV/c 2 ) mass window around the 7r° peak. The 
7T° daughter photons are required to have energies greater than 50 MeV. The 7r°s are then 
constrained to the nominal tt° mass [[J. Candidate rj mesons are required to have invariant 
masses within ±2.5<r of the 77 peak, where a is 10.6 MeV/c 2 and 3.4 MeV/c 2 for the 77 and 
7T + 7r~7r° modes, respectively. For the 7r + 7r~7r° mode, the h + tt~ pair is constrained to a vertex. 
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Both photons from the r] — > 77 mode are required to have E 1 > 100 MeV and the energy 
asymmetry of the daughter photons, ^ 71 ^ 72 ^ , is required to be less than 0.8. We remove ri 

- C/ 71+- C/ 72 

candidates if either of the daughter photons can be combined with any other photon with 
Ej > 100 MeV to form a 7r° candidate. The rj candidates are further constrained to the 
known 77 mass H. Candidate uj mesons are constructed from 7r + 7r~7r° combinations where 
the 7r + 7r~ pair must form a vertex, and the CM momentum of the 7T° is required to be greater 
than 350 MeV/c to reduce the large combinatorial background from low energy photons. 
The invariant mass of the 7r + 7r~7r° combination is required to be within ±30 MeV/c 2 of the 
nominal uj mass || (the natural width of the uj is 8.9 MeV/c 2 ). 

We reconstruct D° mesons in the D° — > K~7T + , K~7t + tt , and i^~7r + 7r~7r + decay modes. 
The CM momentum of the 7r° from D° — ► K^-k + -k° decay is required to be greater than 
300 MeV/c. The invariant mass of the D° candidate is required to be within ±2.5<r of the 
measured D° mass where a, the D° mass resolution, varies between 5.5 and 13 MeV/c 2 de- 
pending on the decay mode. A mass and vertex-constrained kinematic fit is then performed 
on the D° candidates. D*° candidates are reconstructed in the D*° — > D°tt decay mode; 
the minimum photon energy requirement is reduced to 20 MeV for this case. The mass 
difference, 5m = M(D°n°) - M(D°), is required to be within ±2.5cr {a = 0.82 MeV/c 2 ) of 
its nominal value. 

We combine D°s or D*°s with h° meson candidates to form B° candidates. Two 
kinematic variables are used to identify signal candidates, the beam-constrained mass 



M bc = yJ(E^J 2 - (pg M ) 2 and the energy difference AE = E™ - E™ m , where E™ and 
p§ M are the CM energy and momentum of the B° candidate, and E™ m = = 5.29 GeV. 

The typical M& c resolution is 3 MeV/c 2 ; the AE resolution ranges from 17 to 35 MeV, de- 
pending on the decay mode. When more than one -B candidate is found in an event, the 
candidate with the minimum \ 2 is chosen, where x 2 — xjjo + Xho(+Xs m )- Here xjjo is from 
the kinematic fit to the D°, x\o is from the kinematic fit to the h° for the ir° or the rj, while 
Xho = (A(M U! ) 1 'cr(M u; )) 2 for the uj. Here A(M £J ) is the mass difference between measured 
and nominal mass values and a(M u ) is the measured resolution. For the D*°h° modes, ximi 
defined as (A(6m) / a(5m)) 2 , is included in the best candidate selection. 

The background from continuum e + e~ — > qq production is suppressed in the following 
way. We form a Fisher discriminant |J containing seven variables that quantify event 
topology. The Fisher variables include the angle between the thrust axis [|K| of the B 



candidate and the thrust axis of the rest of the event (cos^r), the S±_ variable ||11|| , and five 
modified [[EJ Fox- Wolfram moments |13[| . In the D*°tt and D°u> modes, helicity provides 
additional discrimination. For the D*°7T mode, we define W D .o as the cosine of the angle 
between the D*° flight direction and the direction of the 7r° in the D*° rest frame and require 
|Wd*o| > 0.4. For the D°u mode, we define Ti^ as the cosine of the angle between the B 
flight direction and the normal to the u decay plane in the uj rest frame. We also define a 
variable A as the absolute value of the cross product of the two charged pion momentum 
vectors, \P n + x P n -\, in the uj rest frame. To suppress the qq background, we take the B 
flight direction and the Fisher discriminant, and for modes with uj mesons we include the 
variable A. For the D°uj mode, we also include the variable H w . These variables are then 
combined to form signal (S) and background (BG) probability density functions (PDFs). 
Signal PDFs are determined using Monte Carlo (MC), and background PDFs are obtained 
from Mt, c sideband data. The PDFs are multiplied to form a signal (background) likelihood 
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£s(bg)> and a selection is applied on the likelihood ratio £s/(£s + £bg)- This requirement, 
which has a typical efficiency of 70%, removes more than 90% of the qq background. 

In addition to the qq background, there are large background contributions from color- 
favored B — > DW(rar) - decays and feed-down from D*°h° to D°h° modes. The color-favored 
B — > D(*\nn)~ events give rise to two kinds of backgrounds: those with the same final state 
particles as signal events, and those with missing or extra particles. The former mimics the 
signal distributions in both Mf, c and AE, and selection cuts on variables other than these 
(discussed below) are needed to suppress these backgrounds. The backgrounds with different 
final state particles arise when a pion is missed or a spurious pion is added. In this case, 
the background events may not be distinguishable in the M& c distribution if the missing or 
extra pion has very low momentum, but the AE distribution provides a useful discriminant 
because of the missing pion rest-mass energy. 

The B° — > D* + p~ mode can give the same final states as D°lj and D°rj. However, more 
than 99% of the background events are removed by the uj and r\ mass cuts since the invariant 
mass of the three pions rarely falls within the to or rj mass region. 

The B- -> L>M V final state can contaminate the D^ Q n° mode if the p decay produces 
a fast 7T°. The background events from D°p~ to D°7r° (D*°p~ to D*°7r°) have AE values 
that are below the signal region because of the slow 7r~ that is missed in the reconstruction. 
However, the signal contamination is not negligible due to the large branching fraction of 
B- -> D^°p-. The D°p- mode can also cause a background to D*°rr if a low momentum 
7T° is used to form a D*° candidate. This case is similar to the case where the final state 
particles are the same, and causes a few events in the signal region. About half of these 
events are removed by rejecting events that can be reconstructed as B~ — > D^*'°p~ . This 
reduces the systematic error when fitting for the signal yields and only removes a few percent 
of signal events. 

The B- -> D^°p~ mode can also contaminate the D^ ij channel if a photon from 
the fast 7r° is combined with another photon to form an rj candidate. The contributions of 
these backgrounds in the rj channel, as well as the feed-across from the D^°ti mode, are 
minimized by the 7r° veto described earlier. 

We also check the background contributions from B~ — > D^*'°p'~ (p'~ — > ujtx~) decays 



that have been reported recently by the CLEO Collaboration ||14j| . This two-body decay 
produces high momentum D^°s and us that may fake signal events. Monte Carlo studies 
indicate that the remaining background events are shifted in AE by more than the mass 
of the missing pion and thus are distinguished from signal events by the fit to the AE 
distribution. The D°p'~ mode can also contaminate the D*°u mode if the n~ from p'~ is 
replaced by a 7r°. However, since the 7r~ from p'~ decay carries sizable momentum, the 
kinematics of the final state particles are different from that of the D*°u> signal, and the 
expected background is small. 

The AE distributions for the various D^*'°h° decays, after the application of all selection 
requirements and with M^ c between 5.27 GeV/c 2 and 5.29 GeV/c 2 , are shown in Fig. |2|. The 
plots are fitted to the signal and background functions using a binned maximum likelihood. 
The signal shape is an empirically determined parameterization [15] with parameters ob- 
tained via MC. We observe good agreement between data and MC for the AE distributions 
of color-favored decays such as D* + p~ and D°p~ . The background functions include a com- 
binatorial component and a color-favored component. For the D°h° modes, we also include a 
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component for feed-down from color-suppressed D*°h° modes. The contribution from D°h° 
to D*°h° modes is found to be negligible. The combinatorial component is taken to be 
a first-order polynomial with a slope determined from the AE shape of the M^ c sideband 
(5.20 GeV/c 2 < M bc < 5.26 GeV/c 2 ) data. The shapes of the color-favored and feed-down 
components are modeled by MC. The area of the feed-down component from D*°h° to D°h° 
is fixed by the obtained signal yield from D*°h° with the estimated feed-down efficiency in 
MC. For the D^*'°u mode, the D^*'°p'~ components are also fixed by MC using the measured 
branching fractions |I4 ] . The normalizations of the signal and background components are 
free parameters. 

Table I lists the signal yield, statistical significance, reconstruction efficiency, and branch- 
ing fraction for each D^h mode. In addition, the estimated number of events in the region 
of -0.1 GeV< AE < 0.1 GeV due to backgrounds from generic BB decays, D*°h° feed-down, 
and qq production are given. The systematic errors due to fitting are obtained by varying 
the parameters of the fitting functions within la of their nominal values. The change in the 
signal yield from each variation is added in quadrature to obtain systematic errors from the 
fit. These are between 5% and 15% depending on the decay mode. The statistical signif- 



icance is defined as y — 21n(£(0)/£ max ) where £ max is the likelihood at the nominal signal 
yield and C(0) is the likelihood with the signal yield fixed to zero. We observe signals for 
B — > D°n°, D*°tt°, D°T], and D°uj decays with more than 4cx significance. Independent fits 
to the distributions, after subtracting BB background, confirm these results. We find 
evidence for B° — ► D*°rj and D*°uo with more than 3a significance. For decay modes with 
significance less than 4a, we give 90% confidence level upper limits (UL) on the signal yields 

(iV s UL ) from the relation C(N S ) dN s / f™C(N s )dN s = 0.9, where C(N S ) denotes the 
maximum likelihood with the signal yield fixed at N§. 

The efficiencies for each mode are obtained by Monte Carlo and calibrated by a detailed 
study of detector performance. In particular, studies of tracking, 7r° detection, and parti- 
cle identification give systematic errors in the detection efficiencies of low momentum tt°, 
energetic n°, 77, and to mesons of 10.7%, 7.3%, 9.6%, and 9.5% respectively. The system- 
atic error on the D° meson reconstruction efficiency is estimated to be 12.7% by comparing 
the observed yield of B — > D°7T~ events with the expected yield using known branching 
fractions |5J. A 5% systematic error on the likelihood ratio requirement is determined by 
applying the same procedure to the B~ — > D* tt~ sample and comparing the effects on data 
and MC. The final systematic errors of the branching fractions include the errors in fitting, 
reconstruction efficiency, cut efficiency for background suppression, and the number of BB 
pairs. Assuming the number of B°B° and B + B~ pairs to be equal, we calculate the branch- 
ing fractions for various decay modes given in Table I. The branching fraction upper limits 
are calculated by increasing Ng L and reducing the efficiency by their systematic errors. 

Our measurement of B(B° — ► D°7t°) = (3.1 ± 0.4 ± 0.5) x 10~ 4 is above the previous UL 
of 1.2 x 10~ 4 obtained by CLEO [|TJ. It is also considerably higher than the factorization 
prediction of 0.7 x 10~ 4 ||. This could indicate the presence of final state interactions, or 
other corrections to factorization. It is customary to decompose the B~ — > D°tt~ and 
B° — > D + n~, D°ir° decay amplitudes into isospin | and | components. Our measurement 
of B(B° — > D°tt°), together with the known branching fractions of the other two modes, 
suggest a rescattering phase difference between isospin | and | amplitudes that is 31° ± 7°. 
A similar value, 32° ± 8°, is obtained for the B — > D*tt system. 
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Our measurements of B(B° -> D°i]) = (1.4 ±g;t ± 0.3) x 1(T 4 and -6(5° -»■ D°c;) = 
(1.8 ± 0.5 to' 4 ) x 10~ 4 are also higher than the factorization predictions of 0.5 x 10~ 4 and 
0.7x 10 -4 ||. The central values of the two less significant modes show the same pattern. The 
results for these modes cannot be accommodated by the aforementioned elastic rescattering 
phase. 

In summary, using 23.1 million BB events collected with the Belle detector, we report 
the first observations of color-suppressed B° — > D°tt°, D* tt°, D°rj, and D°u decays, and 
evidence for B° — > D*°r] and D*°uj modes. All these modes have similar branching fractions 
with central values between 1.4 xlCT 4 and 3.1 xl0~ 4 , as given in Table I. They are all 
consistently higher than recent theoretical predictions based on the factorization hypothesis. 
This may be accounted for by additional corrections to the factorization models, or by non- 
factorizable effects such as final state interactions. 
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Resources; the Department of Science and Technology of India; the BK21 program of the 
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TABLE I. The obtained signal yields, statistical significances (S), background events from 
generic BB decays {BB bg), D*°h° feed-down (D*°h°), and qq production (qq bg), reconstruction 
efficiencies (e) (including the sub-decay branching fractions), branching fractions (£>), and 90% 
confidence level upper limits (UL) for B° — > D*°r] and are listed together with the theoretical 

predictions (Th) [2]. The background events are estimated under the signal peak (in the region of 
-0.1 GeV< AE < 0.1 GeV). 

Mode Signal Yield S BB bg D*°h° qq bg e(%) B (x 10~ 4 ) UL (xKT 4 ) Th (x!0~ 4 ) 



L>°7T 126.2 1^1^ 9.3 26?7 L3 145.6 1.79 3.1 ±0.4 ±0.5 0?7 

Z?*°7r° 26.4 4- 1 5 - 9 " 10 - 4 °- 42 2 - 7 la? -ai " L0 

L>V 22.1 1^ ±^ 42 34 07 19.1 0.67 1.4 ± 0.3 1 (h5 

L>*°?7 7.8 ±|;g ± 0.7 3.3 1.4 - 1.5 0.17 2.0 ±g;| ± 0.4 4.6 0.6 

DK> 32.5 tit It" 4.4 5.3(2.3)* LI 58.5 0.80 1.8 ± 0.5 - 07 

D*°uj 16.1 1^ ±2.4 3.0 5.3(1. 5) f 13.8 0.23 3.1 ±\% ± 0.8 7.9 1.7 



* For decay modes with r\ mesons, about 70% of the signal is reconstructed in the 77 channel and the 
remainder is in the 7r + 7r~7r° channel. ' The values in parentheses are the estimated background from 
B~ -» D^°p'- decays. 




d * * d d < < d 

FIG. 1. The external (left) and internal (right) spectator diagrams for B — > Dir decays. 
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AE (GeV) 

FIG. 2. The AE distributions for (a) D°tt°, (b) D*°ir°, (c) D°7], (d) D*°7], (e) D°u, and (f) 
D*°uj. The solid lines show the fit results, the dotted lines are the signals, and the dashed lines show 
the sum of the feed-across and the combinatorial background, with the latter shown separately as 
the dash-dotted lines. 
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